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We propose two dimensional x-ray coherent correlation spectroscopy (2DXCS) for the study of 
interactions between core-electron and valence transitions. This technique might find experimental 
applications in the future when very high intensity x-ray sources become available. Spectra obtained 
by varying two delay periods between pulses show off-diagonal cross-peaks induced by coupling of 
core transitions of two different types. Calculations of the Nls and Ols signals of aminophenol 
isomers illustrate how novel information about many-body effects in electronic structure and exci- 
tations of molecules can be extracted from these spectra. 

PACS numbers: 33.20.Rm, 42.65.Re 



X-ray absorption spectroscopy (XANES, EXAFS)[T] 
and its time-resolved extensions [3J 2] provide a di- 
rect probe for electronic structure of molecules with sub- 
atomic and subfemtosecond resolution. Novel ultrabright 
x-ray sources such as future free-electron laser (XFEL) 
[5] or high-harmonic (HH) sources [B] may make it possi- 
ble to perform nonlinear experiments with multiple x-ray 
pulses. All-x-ray nonlinear signals could provide more de- 
tailed information on the electronic structure and dynam- 
ics than available from time-resolved XANES, by prob- 
ing states with multiple core electrons excited. Many 
proposed applications of the new sources make use of 
their ultrashort temporal resolution and high intensity to 
monitor dynamical processes in real time. Two-photon 
absorption [7 and x-ray driven molecular dynamics |5] 
have been demonstrated using HH sources. Techniques 
such as diffraction or pump probe do not rely on the 
coherence properties of the beams. The technique con- 
sidered in this Letter, in contrast, depends also on pulse 
coherence in an essential way, and should become feasible 
once high intensity, attosecond |HJ E01 EH dj transform- 
limited pulses become available. Such pulses should allow 
to control and manipulate the coherence of core excita- 
tions and use it as a window into correlations between 
different regions of the molecule. Similar ideas are effec- 
tively used in multidimensional NMR spectroscopy [13 
to probe correlations between spin dynamics in controlled 
time periods using elaborate pulse sequences. The sig- 
nals are interpreted in terms of multiple correlation func- 
tions which provide fundamentally new types of informa- 
tion compared to one dimensional techniques. The same 
ideas were recently extended to the infrared and optical 
regimes [H ESI EES HZ] ■ 

In this letter we propose a new class of two dimensional 
x-ray coherent correlation spectroscopy (2DXCS) tech- 
niques and demonstrate how they could provide a unique 
probe for interactions between the core-transitions and 
electronic states that mediate these interactions. In- 
frared femtosecond 2D techniques can excite molecular 
vibrations impulsively and probe the subsequent corre- 



lated dynamics of nuclear wavepackets. Similarly at- 
tosecond x-ray pulses resonant with core transitions can 
excite valence electrons impulsively and probe correla- 
tions in dynamical events of resulting electron wavepack- 
ets. Since core transitions are highly localized to the ab- 
sorbing atoms, these techniques also offer a high spatial 
resolution. 

We consider a time-resolved coherent all-x-ray four- 
wave mixing process carried out by subjecting the 
molecule to a sequence of three pulses. The first pulse 
has wavevector ki and carrier frequency oj\, followed se- 
quentially by the other two pulses k^,uj2 and &3, lu^. 
The signal is heterodyne-detected with a fourth pulse 
^4,^4. The coherent nonlinear response generated in the 
kj = k\ + ki + &3 phase matching direction is recorded 
as a function of the delays t\ , t 2 , and t 3 between consec- 
utive pulses. 2DXCS is obtained by a Fourier transform 
of the signal with respect to two delays, and displaying 
it as a two-dimensional frequency correlation plot 

S! 3) (ft 3 ,t 2 ,tti) = JJ^dhdhS^ihMMV^e 1 ^. 

Like in NMR, transitions of the same type contribute 
to the diagonal part Q\ = fl 3 of the correlation spectrum, 
while features arising from interactions among spectrally 
separated transitions appear as off-diagonal cross peaks. 
By spreading the signal over multiple frequency axes, the 
weak signatures of interactions between different transi- 
tions can be separated from strong same-transition sig- 
nals [HE!]- 

We illustrate the power of these techniques for Nls and 
Ols transitions of aminophenols (Figure [TJ. The polar- 
ization induced in the molecule by interactions with the 
x-ray pulses is calculated in terms of nonlinear response 
functions. Different time orderings appearing in the per- 
turbative expansion of the density matrix are described 
as Liouville space-pathways and represented by double- 
side Feynman diagrams 19J. Assuming temporally well 
separated pulses, the signal within the rotating-wave ap- 
proximation (RWA) where we only retain the resonant 
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FIG. 1: (a) Para and ortho isomers of aminophenol. (b) 
Valence and core-excited states of aminophenols: g denotes 
states with no core electron excited (including the ground 
state go), ejv (eo) denotes states with the Nls (Ols) electron 
excited, and / denotes states with both Nls and Ols electrons 
excited, (c) Double-side Feynman diagrams representing the 
two contributions to the cross peak 2DXCS signal [Eq. [I] 



contributions is given by 
S r 3 (ta, t 2 , ti) 



(1) 



(2) 



where the response function 

Rf\h,t 2 ,h) 

=i 3 [(Bi(0)B+(t 2 + h)B±{t 3 + h + t^B+ih)) 

+ (%;(0)j%(ti)Bz(t 3 + t 2 + h)B+(t 2 + h)) 

- (Bf (0)B^(t 3 + t 2 + h)B+(t 2 + t^B+ih))] 



is expressed in terms of the core exciton creation and 
annihilation operators weighted by the pulse envelopes 
and expanded in terms of the molecular electronic states 



Here, £j(u>) is the Fourier transform of the complex en- 
velope of the jth pulse £j{uj) = J dtEj{r - r^e'^-^ . 
and ujj its carrier frequency. 

In the following simulations we assumed the pulse 
bandwidths (10 eV for a 125 attosecond pulse) to be 
much smaller than the splitting between the Nls and 
Ols transitions (~ 120 eV). The carrier frequencies uij 



thus select the desired type of core transition, while the 
pulse envelopes £j(ui) control which valence excitations 
of the selected type contribute to the response within the 
pulse bandwidths. 

The diagonal part of the 2DXCS signal is obtained by 
tuning all four uij either to Nls or to Ols core transi- 
tions, while the cross peaks are obtained in the sequen- 
tial two-color pulse configuration where u>\ — oj 2 = ujn 
are tuned to the Nls and co 3 = U4 = uoo to the Ols 
transition. In the latter case, the signal consists of the 
ground-state bleaching (GSB) and excited-state absorp- 
tion (ESA) terms shown in Figure Id. The former repre- 
sents the effect of molecules missing in the ground-state 
after two interactions with lun and cannot absorb the too 
photon. The latter represents the additional absorption 
of molecules from the singly- to doubly-excited states. 
The signal is given by 
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where ^, a b and T a t, are the transition dipole and the de- 
phasing rate of the ab transitions, respectively. 

The cross peaks are intrinsically nonadditive and 
highly sensitive to the coupling between core transitions. 
When the two core transitions are decoupled, the cross 
peaks vanish identically due to the interference between 
the one (GSB) and two-core exciton (ESA) terms. If 
the core transitions are independent, we can express the 
two-core exciton states as direct products |/) = |eie 3 ), 
u>e 3 g — LjJ fe 1 , and the two denominators in Eq. Q become 
identical. Similarly the transition dipole elements will be 
the same and the two terms which have an opposite sign 
exactly cancel. Having distinct features in the spectra 
that are induced by correlations is the main merit of 2D 
techniques. In NMR, cross peaks are used to extract spin 
coupling and convert this information into the structure 
of complex molecules. 

First-principles simulations of the cross peaks depen- 
dence on the states with two core electrons excited are 
particularly challenging. We describe the core transition 
as a response of N+l valence electrons to an instantly 
switched core hole using the Nozieres - De Dominicus 
Hamiltonian [20] . The core-hole potential was taken into 
account by incrementing the charge of the corresponding 
nuclei (equivalent-core or Z+l approximation |21jV The 
valence excitations of the equivalent-core molecule were 
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FIG. 2: (Color online) Simulated Ols/Nls 2DXCS cross peak (1,2 = 0) and its GSB and ESA components for para- (upper 
panel) and ortho-aminophenol (lower panel). Insets show the simulated Ols and Nls XANES plotted as a function of u — Uj. 



represented by unoccupied Kohn-Sham orbitals. This al- 
lowed us to maintain a simple correspondence between 
the molecular orbital picture and the many-body states, 
which facilitated the connection between their properties 
and the 2DXCS spectra. 

Figure [2] shows the total Ols/Nls cross peak 

S { I ~ 3 \ni,t 2 = 0,fi 3 ) [Eq. (0] and its GSB and ESA 
components for para and ortho isomers of aminophcnol. 
The pulse carrier frequencies oj\ and L02 were tuned to 
the lowest Ols transition (^535 eV), W3 and 0J4 to the 
lowest Nls transition (~401 eV). We assume the 10 eV 
pulse bandwidths of (±5 eV around the carrier frequen- 
cies). The orbitals of the original and equivalent-core 
molecules were calculated using the B3LYP exchange- 
correlation functional and the 6-311G** set of Gaussian- 
type atomic orbitals. 

In both isomers, the calculated GSB component is 
identical to a two-dimensional product of the Ols and 
Nls XANES spectra. In the para isomer, the O and 
N atoms are spatially separted, and promoting the Ols 
electron only weakly affects the Nls transitions. Conse- 
quently, the ESA and GSB components are similar and 
the total 2DXCS cross peak is weak. In the orho isomer, 
the O and N atoms are close, and exciting the Ols elec- 
tron strongly affects the Nls transitions. The ESA peaks 
are shifted, resulting in a stronger 2DXCS cross peak. 
Thus, the 2DXCS cross peak is highly sensitive to the 
relative position of the N and O atoms in the aminophe- 
nol isomers. 

Figures [3] shows the features of the simulated Ols 



XANES and Ols/Nls 2DXCS cross peak of para- and 
ortho-aminophenol in the region marked on Figure [2] 
The intensity scale was magnified by a factor of two to 
better show the structure of each peak. The equivalent- 
core molecular orbitals responsible for the three stronger 
peaks (marked A, B, and C) in the Ols XANES are 
shown as well. For example, orbital A is populated by 
the promoted Ols electron in the lowest Ols transition 
and, in the single-orbital approximation, represents the 
lowest Ols core-excited state of para-aminophenol. 

The XANES signals are insensitive to the precise form 
of the wavefunctions of the corresponding states. For ex- 
ample, despite the difference in the shape of orbital B 
between two isomers, its profile in the vicinity of the O 
atom remains the same and so is the corresponding peak 
in XANES. However, its contribution to the 2DXCS cross 
peak is markedly different. In the para isomer, orbital B 
is delocalized and extends over both O and N atoms. 
Correspondingly, promoting the Ols electron to this or- 
bital significantly affects the Nls transition, resulting in 
a distinguishable cross peak. In the ortho isomer, orbital 
B vanishes in the vicinity of the N atom. Promoting the 
Ols electron does not affect the Nls transitions and the 
cross peak vanishes. The cross peak thus carries informa- 
tion about the wavefunctions of these states. Similarly, 
despite a stronger dipole coupling, the contribution of or- 
bital A of para-aminophenol to the cross peak is much 
weaker than that of orbital B, due to the localization of 
the former to the O atom. In the ortho isomer, however, 
the corresponding cross peak is strong, indicating that 
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FIG. 3: (Color online) Simulated Ols XANES and Ols/Nls 
2DXCS cross peak (at ti — 0) of para- (left panel) and orhto- 
aminophenol (right panel). Insets show the molecular orbitals 
populated by the promoted Ols electron for each of core- 
excited states contributing to the signal. 



the atoms are sufficiently close for the Nls transitions to 
be affected even though the Ols electron remains in the 
vicinity of the O atom. 

In conclusion, we have demonstrated how the off- 
diagonal part of the 2DXCS spectra depends on the in- 
teractions between core transitions, providing informa- 
tion about the electronic states that mediate this in- 
teractions. Simulations of the cross peak 2DXCS sig- 
nal of aminophcnol showed that it provides further in- 
sights into the wavefunctions of these states beyond the 
dipole coupling available from XANES. The dependence 
of the 2DXCS spectra on transitions of different char- 
acter (valence, core transitions with one and two core 
electrons excited) makes it particularly challenging for 
first-principles simulations. Predicting and interpreting 
the 2DXCS spectra will provide a critical test for elec- 
tron structure computations and should stimulate the 
development of accurate methods for excited states of 
two-electron character. 

The 2DXCS measurements proposed here are beyond 
the capabilities of existing x-ray sources. Generating 
the required pulse sequences with sub femtosecond du- 
ration, high flux and controlled timing and phase will re- 
quire to overcome many experimental difficulties. Given 
the calculated values of the N Is and O Is transition 
dipole moments of 0.1 atomic units, the laser focal area 
of 1CP 10 cm 2 , and the surface density on the order of 10 14 
aminophcnol molecules per cm 2 , we estimate the absorp- 
tion probability to be on the order of 10 -5 . In the fully 
resonant four-wave-mixing signal there is thus 1 gener- 
ated photon per 10 15 incoming photons. Even with the 
XFEL capability of generating 10 13 photons per pulse, 
high repetition rates and long acquisition times will be 



necessary to make the experiment feasible. Other experi- 
mental challenges include reducing the duration of XFEL 
pulse [12 , mainting longitudinal coherence, and gener- 
ating multiple pulses with controlled timing. The HH 
sources have the required pulse duration and are fully 
spatially coherent, but have much lower brilliance [BJ. 
New phase matching technique that works at high pho- 
ton energies are being developed [22] . This paper aims at 
establishing the theoretical basis only, and as it took 20 
years to extend this NMR radiowave technology to the 
infrared and the visible regions, similar time and break- 
throughs will be needed to realize these measurements in 
the x-ray regime. 
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